Thermodynamic and transport properties of NiCl 2 -4SC(NH 2 )2: role of strong mass renormalization 
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Several quantum paramagnets exhibit magnetic field-induced quantum phase transitions to an antiferromag- 
netic state that exists for H c \ < H < H C 2. For some of these compounds, there is a significant asymmetry 
between the low- and high-field transitions. We present specific heat and thermal conductivity measurements in 
NiCl2-4SC(NH2)2, together with calculations which show that the asymmetry is caused by a strong mass renor- 
malization due to quantum fluctuations for H < H c i that are absent for H > H C 2- We argue that the enigmatic 
lack of asymmetry in thermal conductivity is due to a concomitant renormalization of the impurity scattering. 



PACS numbers: 75.10.Jm, 75.40-s, 75.40.Cx 

The correspondence between a spin system and a gas of 
bosons has been very fruitful for describing field-induced or- 
dered phases in a large class of quantum paramagnets 03451 . 
In this analogy, a magnetic field H plays the role of the chem- 
ical potential, which, upon reaching a critical value H c \, in- 
duces a T = Bose-Einstein condensation (BEC), provided 
that the number of bosons is conserved, the kinetic energy is 
dominant, and the spatial dimension d > 1. Such a BEC state 
corresponds to a canted XY magnetic ordering of the spins. 

At the BEC quantum critical point (QCP), the low-energy 
bosonic excitations have a quadratic dispersion u> = k 2 /2m*, 
where m* is the effective mass. This mass is renormalized by 
quantum fluctuations in the paramagnetic phase H<H C % . In 
magnets with H c i <C H C 2 the renormalization can be expected 
to be very strong because of the proximity to the magnetic 
instability. The transition at H c % should be contrasted with 
the second BEC-QCP that takes place at the saturation field 
H C 2 (6J. Since the field induced magnetization is a conserved 
quantity, there are no quantum fluctuations and no mass renor- 
malization for the fully polarized phase above H C 2, i.e., the 
bare mass m can be obtained from the single-particle excita- 
tion spectrum at H > H C 2- Thus, quantum paramagnets are 
ideal for studying mass renormalization effects because the 
effective and the bare bosonic masses can be obtained from 
two different QCP 's that occur in the same material. 

Here we present theoretical and experimental evidence for 
a strong mass renormalization effect, m/m* ~ 3, in MQ2- 
4SC(NH 2 ) 2 [referred to as DTN]. We will show that the large 
asymmetry between the peaks in the low-temperature spe- 
cific heat, C V (H), in the vicinity of H c \ and H C 2 is closely 
described by analytical and Quantum Monte Carlo (QMC) 
calculations. The mass renormalization also explains simi- 
lar asymmetries observed in other properties of DTN, such as 
magnetization [7 1, electron spin resonance [8 1, sound velocity 



||9l [Toll , and magnetostriction [fT|. In a remarkable contrast 
to these properties, peaks in the low-temperature thermal con- 
ductivity, k, near H c \ and H C 2 do not show any substantial 
asymmetry. We provide an explanation to this dichotomy by 
demonstrating that the leading boson-impurity scattering am- 
plitude is also renormalized by quantum fluctuations, effec- 
tively canceling mass renormalization effect in k. 

DTN is a quantum magnet with tetragonal crystal symme- 
try that exhibits a field-induced BEC iTTl l8l \T2MU\ to a very 
good approximation |15|. The dominant single-ion uniaxial 
anisotropy D = 8.9K splits the Ni 5 = 1 triplet into an S z = 
ground state and an S z — ±1 excited doublet. The antiferro- 
magnetic exchange coupling between Ni ions is J c = 2.2 K 
along the c-axis and J a = 0.18 K along the a- and b-axes, 
while the gyromagnetic factor along the c-axis is g = 2.26 
|8l . A magnetic field applied along the c-axis lowers the 
energy of the S z = 1 state producing a T = BEC transi- 
tion at H c i =2.1 T. The long-range order occurs in a dome- 
shaped region of the T — H phase diagram between H c \ and 
H C 2 = 12.5 T and below the maximum ordering temperature 
T max ~ 1.2K lHH. The T 3/2 dependence of the critical field 
expected for a BEC-QCP has been established via direct mea- 
surements of the phase boundary with ac susceptibility down 
to 1 mK 1 14 1, and by magnetization measurements lfT6ll . The 
asymmetry between H c \ and H C 2 I17 UT21 can also be seen di- 
rectly in the skewed shape of the phase diagram [ 12 1. 

The Hamiltonian describing the 5 = 1 spin degrees of free- 
dom of DTN in external field is given by (8) 

7i = J - S r • S r+e„ + D £ (S Z f - #r , (D 

r ,v r r 

where e„ are the primitive vectors of the lattice, v = {a, b, c}, 
and h = g/iBH. We introduce Schwinger bosons associated 
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with the fundamental representation of SU(3) that obey the 
constraint Yl m b rm b rm — 1. The subscript m = {I, 0,t} 
labels the eigenstates of with the eigenvalues { — 1, 0, 1}. 
The spin operators in this representation are: 



lbt t b r0 + bt b ri )A2) 



with n rm = b\ m b Tm . We enforce the constraint by introduc- 
ing spatially uniform Lagrange multiplier /i 



H = H + /i ( b U b U + b tA + b lo b to 1 



(3) 



The lowest energy state in the H < H cl paramagnetic regime 
is &r |0) and the ground state corresponds to a non-zero ex- 
pectation value of the S z = boson: b\ = b r0 = s. By using 
the spin representation |2} with the mean-field value for b^ 

and neglecting higher-order terms in powers of &£m> we °b~ 
tain the Hamiltonian in the harmonic approximation 



n = E Q 
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with Ak a = (/i + s 2 ek — h a ) and £>k = s 2 ek, where E — 
N(fi—D)(s 2 —1) is the bare ground-state energy, N is the num- 
ber of sites, a = {t, I}, h a = ±h, <j = —o, S^J are the Fourier 
transformed bosonic operators, and 6k = 2 ^ Jv cos fc„. The 
anomalous terms indicate that bosons with opposite S z are 
created and annihilated in the ground state. These are the 
quantum fluctuations that lead to renormalization of the quasi- 
particle dispersion relation. The Hamiltonian Q is diagonal- 
ized by the Bogolyubov transformation 



b ka = u k /3 k(T + w k /?! 



(5) 



where w k u k = 13^/2^, 



v i = (/•* + s2e k) and 



^ y? + 2^s 2 e k - The resultant diagonal form of H is 



n = E () 



k 



(6) 



Thus, the low-energy spectrum for h < h cl is uj^ = w£ — h. 
The band uj^ has a minimum at the antiferromagnetic wave- 
vector Q = (tt, it, 7r) with the gap A< = cjq — h, whose van- 
ishing point defines the critical field h c \ = g^sHd = Wq. The 
ground state energy is also affected by quantum fluctuations: 



En — En 



(7) 



The saddle point conditions, dEo /ds = dEn /dpi = 0, lead to 
the self-consistent equations for the parameters s and /i, 



s z = 2- 



k k 



, D = n 



N 
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(8) 



Using the Hamiltonian parameters for DTN given in Ref. |8], 
the resulting values are s 2 — 0.92 and // = 10. 3K. 



This low-energy theory is valid only for H < H c ±. For 
H > H C 2 spins are fully polarized and the spectrum can be 
computed exactly. Since there are no quantum fluctuations 
for H > H c2 , the exact value of h C 2 is h C 2 — QUeHcI = 
D — 2eq, while the unrenormalized excitation spectrum is 
w^* = e k — eQ + h — h C 2, which also has a minimum at Q with 
the gap A > = h — h C 2- Since only the excitations near k = Q 
are important at low temperatures, we define the mass tensors 
for H < H c \ and H > H c2 as: 
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Then the mass renormalization factor is given by 
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(9) 



(10) 



For the parameters of the Hamiltonian from Ref. [8 1, we ob- 
tain m vv jm* vl , ~ 3.2. Such a large difference of masses must 
readily demonstrate itself in the strong asymmetry of the C v 
vs H curves near H c i and H C 2 as well as in the slopes of 
the specific heat dependence on T at the critical fields where 
C v oc (Tm) 3 ' 2 . These theoretical expectations are supported 
by the experimental C V (T, H) data shown in Figs. [T] and [2] 

The C p was measured in single crystals of DTN grown 
from aqueous solutions of thiourea and nickel chloride, with 
magnetic field applied along the crystalline c-axis. The ex- 
perimental C p vs H was obtained using an AC technique 
ifTTIl . while sweeping the magnetic field in a 3 He fridge fur- 
bished with a 17 T superconducting magnet system at the Na- 
tional High Magnetic Field Laboratory (NHMFL) and the Los 
Alamos National Laboratory. We also used the standard ther- 
mal relaxation method to obtain C v vs T with a dilution re- 
frigerator in a 16 T Physical Properties Measurement System 
at Quantum Design, Inc. A strongly asymmetric C p vs H is 
shown in Fig.[T]for fixed temperatures T = 0.75K (green line) 
and T = 0.4K (red line), alongside the results of the QMC 
simulations of the spin Hamiltonian in Eq. ([T]) in a 8 x 8 x 32 
lattice (solid symbols). The agreement between the QMC re- 
sults and the experimental data is very good. For the low- 
est temperature, the asymmetry of the peaks is C P 2/C P \ ~ 6, 
close to (m/m*) 3 / 2 «5.7 expected from the theory above. 

Our C v vs T experimental data close to H c x (top panel) and 
H C 2 (bottom panel) are displayed in Fig. [2] (lines with sym- 
bols). The results of the QMC simulations of Eq.Q at H c i 
and H C 2 are shown by the solid lines. The dashed lines corre- 
spond to the analytical calculation of C v (T) for a dilute gas of 
hard-core bosons [Eq. (|6|] that have the spectra given by 
and at H cl and H c2 . The on-site boson-boson repulsion is 
taken into account at the mean-field level by summing the lad- 
der diagrams (see Ref. [3 1). Since this approach is only valid 
for low density of bosons, it agrees closely with the QMC re- 
sults at low T, but deviates from them at higher temperatures. 
The very good agreement between the theoretical results and 
the experimental curves at H c % and H C 2 confirms quantita- 
tively the expected mass renormalization for H < H c \. 
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FIG. 1 : Specific heat as a function of magnetic field for two temper- 
atures. Measurements and QMC simulations were also performed 
at the other temperatures showing nearly perfect agreement and the 
same characteristic behavior of C p vs H (not included for clarity). 
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FIG. 2: Specific heat data as a function of T for magnetic fields 
near (a) H c \ and (b) H c2 . A Schottky anomaly tail, A/T 2 with 
A — 0.018J/molK 3 , has been subtracted for fields near H C 2- The 
full lines correspond to QMC simulations of Eq.jTJ) for H — H c i(2) 
and the parameters of Ref. (8). The dashed lines are analytical cal- 
culations. 



The thermal conductivity was measured in DTN single 
crystals using the standard uniaxial heat flow method, where 
the temperature difference was produced by a heater attached 
to one end of the sample and monitored with a matched pair 
of RuC>2 thermometers. The heat flow and the magnetic field 
were parallel to the c-axis. Similar observations were reported 
in Ref. 18 although their data at base temperature (380 mK) 
do not agree with ours, measured down to 300 mK. 

The lighter mass of bosons for H < H c \ implies not only 
large asymmetry between the peaks in the specific heat field 



dependence, but also similar asymmetries in a number of other 
properties of DTN flT UTTll . all exhibiting a much stronger 
anomaly at iJ C 2 than at H c %, In contrast, the low-temperature 
thermal conductivity does not show any substantial asymme- 
try between the H c i and H C 2 data. Fig. [3] shows the field 
dependence of the thermal conductivity, k, normalized to the 
H = value, k(0), for several low values of T. Since the 
H = magnetic excitation spectrum has a gap of about 3K 
[8 1, only phonons contribute to «(0) at low temperatures. The 
behavior of k changes qualitatively in the field because the 
gap is closed between H c i and H C 2- The low-temperature 
magnetic excitations provide a substantial contribution to the 
thermal conductivity as is clear from n(H) / k(0) being > 1 in 
Fig. [3] Here we focus on the low-temperature behavior of k 
at the critical points H c i and H c2 . A detailed analysis of the 
other aspects of n will be provided elsewhere |[T9l . 

At low enough temperatures, scattering of bosons on each 
other should diminish because their concentration will be 
small. Consequently, the leading scattering in this regime 
should be due to defects. In the second Born approximation, 
the disorder-averaged inverse mean-free path of an excitation 
of mass m due to scattering on point-like impurities is [20] 



2tt 



\V\ 2 m 2 



(11) 



where rtj is the impurity concentration and V is the effective 
impurity potential. When the excitation gap vanishes at H c x 
or H C 2, thermal conductivity at low-T can be written as 
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K OC 

m 



k 3 dk oc niT 2 ■ £ oc 



T 2 



nim\V\' 2 



(12) 



The theoretical temperature dependence, k oc T 2 , is in a 
good agreement with the measured low-T thermal conductiv- 
ity. However, this dependence does not shed any light on the 
lack of asymmetry between the peaks. Since vacancies or sub- 
stitutional impurities are expected to be rare in clean enough 
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FIG. 3: Thermal conductivity of DTN along the c-axis as a func- 
tion of magnetic field for several low-T values. Inset: Theoretical 
prediction for the peak ratios K2/K1 and C v %jG v \ vs H c \/ H C 2- 
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systems like DTN, we can assume that random lattice distor- 
tions are the most common source of disorder. Because D is 
the largest parameter in Eq.Q, the most significant effect of 
these distortions is a real space modulation of D: 



HP = SD(S?) 2 => SD 



r,k,k' 



e iR.(k-k') 6 t b (13) 



where i is the impurity site and we have used the mapping pi. 

This is where the renormalization due to quantum fluctua- 
tions becomes crucial again. For H = H C 2 the impurity scat- 



tering in ( 13 1 is not renormalized due to the absence of quan- 
tum fluctuations and V2 = SD. On the other hand, for H = H c \ 
the scattering is affected by the quantum fluctuations. Since 
the dressed bosonic excitations are related to the bare ones 
through Eq. d5), this transforms impurity scattering ( fT3) i into 



qjD 
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Thus, the impurity potential at H c i is V\ = 8D (t 



which will modify the mean-free path in (111. After some 
algebra utilizing Eqs. ( fT2) , |5]l and ( [T0| , we finally obtain 



K2 



m ■ £2 _ / m \ 1 / 
m* ■ li Im* / 4s 4 I 




(15) 



This expression contains a large prefactor (m/m*) coming 
from the renormalization of the density of states and veloc- 
ity in ( 12 1, and would formally imply a larger peak at H c2 , 
similar to the specific heat and other quantities. However, 
this effect is partially compensated by the numerical factor 
s» l/4+0((m*/m) 2 ), which comes from the renormalization 
of the mean-free path. By using the DTN parameters, we ob- 
tain Kz/kx s» 1.1 in an excellent agreement with the data in 
Fig. [3] Thus, the mass renormalization effect in thermal con- 
ductivity is compensated by a similar renormalization effect in 
the impurity scattering. To show that this is not a mere coinci- 
dence, we provide our prediction for the H c \ /H C 2 dependence 
of the peak ratios in thermal conductivity K2 / ki and specific 
heat C V 2,/C V \ on H c i/H c2 (see inset in Fig. [5]). Here we used 
the relation between the mass ratio and H c i/H C 2 given by 
Eq. ( 10 1. The vertical line corresponds to the DTN value of 
H c \/H c2 s» 0.17. It is remarkable that k 2 /ki and C^/Cui 
behave in very different ways. In particular, K^/ti ~ 1 while 
the peaks in C v are very asymmetric for 0.1 < H cl /H C 2 < 1. 
With this insight, we also suggest an experimental verification 
of our theory by conducting the heat conductivity measure- 
ment in DTN under pressure. A modest decrease of H c \ by 
IT should lead to an increase in K2/K1 by a factor of 2. 



The leading impurity scattering ( 13 1 and, consequently, the 



resulting expression for the ratio k 2 /ki in (15 1 will remain 
valid for the other BEC magnets even though they may not be 
dominated by the single-ion anisotropy term. For instance, 
in the dimer-based systems |4), the disorder in the leading 
intra-dimer coupling translates into the local modulation of 



the chemical potential which is equivalent to our Eq. (13i 
Thus, our Eq. ( 15 1 can be verified in other BEC compounds. 



In conclusion, by using the example of DTN, we connected 
the asymmetry in the physical properties of BEC magnets with 
the mass renormalization of the elementary excitations due 
to quantum fluctuations of the paramagnetic state. We also 
resolved the enigmatic absence of this asymmetry in the low- 
T thermal conductivity by identifying the leading scattering 
mechanism and by demonstrating that the renormalization of 
the latter compensates the mass renormalization effect. 
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